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Abstract 

The log returns of financial time series are usually modeled by means of the station- 
ary GARCH(1,1) stochastic process or its generalizations which can not properly 
describe the nonstationary deterministic components of the original series. We an- 
alyze the influence of deterministic trends on the GARCH(1,1) parameters using 
Monte Carlo simulations. The statistical ensembles contain numerically generated 
time series composed by GARCH(1,1) noise superposed on deterministic trends. 
The G ARCH (1,1) parameters characteristic for financial time series longer than 
one year are not affected by the detrending errors. We also show that if the ARCH 
coefficient is greater than the GARCH coefficient, then the estimated G ARCH (1,1) 
parameters depend on the number of monotonic parts of the trend and on the ratio 
between the trend and the noise amplitudes. 
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1 Introduction 



The log returns of financial time series {Pt} (share prices, stock indices, foreign 
exchange rates, etc.) 

= log(Pi/Pt-i) = log (Pt) - log (Pi_i) , (1) 
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usually presents the following features: they are uncorrelated, their volatility 
clusters, they have fat-tailed distributions (leptokurtosis), a leverage effect is 
present (changes in stock prices tend to be negatively correlated with changes 
in volatility), their autocorrelation function decays exponentially, their abso- 
lute values present a long range dependence [2]. One of the most used stochas- 
tic models that reproduces some of these features is the Generalized Auto 
Regressive Conditional Heteroskedasticity (GARCH) model having the vari- 
ance expressed as a linear function of past squared innovations and earlier 
calculated conditional variances [T]. There are various generalizations of the 
GARCH model, however the most used in practical applications is its simplest 
form, GARCH(1,1). 

According to relation ([1]) the stationary GARCH(1,1) process is suitable for 
modeling time series for which {log Pt} has a linear trend, i.e. the original time 
series {Pt} contains an exponential trend. But the nonlinear trends in {log Pt} 
are not eliminated by the differentiation ([T]). This problem is amplified in the 
case of a nonmonotonic trend. An alternative to the stationary modeling of 
financial series is the hypothesis of a nonstationary evolution. For example 
Starica and Granger [5] propose a nonstationary model locally approximated 
by a stationary one 

Xt = fi{t) + a {t) St , (2) 

where St are i.i.d with E{et) = and E{ef) = 1 and the unconditional mean 
H (t) and the unconditional variance a (t) are functions of t. If a nonstationary 
series ([2]) is modeled with a stationary process, then the deterministic trend 
fi(t) is confounded with a stochastic trend and the model tends to approach 
its nonstationarity limit. In the case of GARCH model this is the so called 
IGARCH effect. 

In this paper we study the influence of a deterministic trend on the parame- 
ters of GARCH(1,1) model. We use a Monte Carlo method in order to evalu- 
ate the influence of the detrending errors on the variability of the estimated 
GARCH(1,1) parameters. The paper is organized as follows. In the following 
section we shortly present the GARCH(1,1) model and we study the intrin- 
sic variability of its parameters. In the third section an automatic method 
to generate artiflcially trends is described. Then we present the variability of 
the GARCH(1,1) parameters due to the detrending of an artiflcially added 
deterministic component and the last section is dedicated to conclusions. 



2 GARCH(1,1) model 

The GARCH (1,1) process is a real- valued discrete time stochastic process {xt} 
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Fig. 1. Dow Jones Composite index series over the period 01.02.1980-31.12.1999 and 
three index series simulated with G ARCH (1,1) model 

Xt \ {xt-i, crf_i} ~ (O, dj) 
at = K + axl_^ + (3al_^ 

where > 0, a > 0, /3 > [1]. If a + /3 < 1, then GARCH(1,1) pro- 
cess is wide sense stationary with E {xt) = 0, va,T{xt) = K/ {1 — a — [3) and 
cov(xt, Xs) =var(xt) 5ts. 

We analyze the daily Dow Jones Composite (DJC) series containing N = 5089 
values between 01 February 1980 and 31 December 1999 when a deterministic 
trend is likely to exist. The parameters of the GARCH(1,1) model obtained 
using the maximum likelihood method for the log returns of this series are: 
ar^j = 0.0837, Pdj = 0.8898, Kdj = 2.5 ■ 10"^ With these values we have 
generated three index series with the same initial value on the same time 
interval. In Fig. [1] one observes large differences between the generated series 
and the initial one, especially for large t. This behavior is in accordance with 
the fact that GARCH model is suitable only for relatively short periods of 
time |3]. When the series length is large, then the GARCH parameters are 
close to the nonstationarity limit (a-f/3 = 1) and the deterministic trend (if it 
exists) is lost or is strongly distorted because it is replaced with a stochastic 
one. In the case of the analyzed DJC index a^j + Pdj = 0.9735. Therefore, 
the variability of the realizations of a GARCH process with given parameters 
is very large, especially when the series length is large. 

In order to correctly evaluate the variability of the parameters due to detrend- 
ing, it must be compared with the intrinsic variability of the GARCH (1,1) 
parameters for time series without trends. The intrinsic variability is deter- 
mined using a Monte Carlo simulation. For a given length N we generate 
500 realizations of a GARCH(1,1) process with the parameters (a^j, Pdj, 
Kdj). For each realization the GARCH(1,1) parameters are estimated by ap- 
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Fig. 2. The relative standard deviation of the GARCH(1,1) parameters (a) and the 
mean of P (b) for statistical ensembles with 500 numerically generated series with 
DJ parameters for each length A'^. 

plying the maximum likelihood method. The relative standard deviation of 
the GARCH(1,1) parameters decreases to an almost stationary value for large 
values of (Fig. [2^). Generally, the mean of the estimated parameters almost 
coincides with the values used for generating the series with the exception of 
f3 (Fig. In the following we consider series with length = 6000 which 
have a variability near the stationary value. 



3 Automatically generated trends 



In order to be representative, the Monte Carlo statistics must contain a large 
number of numerical simulations with variability comparable with those ap- 
pearing in practical applications. We describe an automatic method for gener- 
ating time series containing a deterministic trend which satisfies these condi- 
tions. The generation of a large number of trends with a significant variability 
using a fixed functional form requires a large number of parameters. For exam- 
ple, a polynomial trend must have a large enough degree, hence the number of 
its coefficients is also large. If we choose the coefficients by means of a random 
algorithm, then the form of the generated trend is difficult to be controlled. 
Usually the resulting trend has only a few parts with significant monotonic 
variation. 

We generate a trend {/„}, n = 1,2,...,N by joining together s monotonic 
semiperiods of sinus with random amplitudes and lengths. In this way we 
obtain a large enough variability for the generated trends and we can control 
the number and the amplitude of its monotonic parts. We need only three 
parameters: the length of the series A^ (in the following A^ = 6000), the number 
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Fig. 3. Artificially generated trends with 2, 3, and 4 monotonic parts 

of monotonic parts s (in our tests s = 1, 2, 3, 4) and the minimum number of 
points in a part equal with 50. The amplitudes of the sinusoidal parts will 
be random numbers with uniform distribution G (0,1). The value of the 
trend at the point n of the part p, Np < n ^ ^p+i, is given by the recurrence 
relation 



fn = fN, + i-iyap 



TT 



sm 



1 + 2 



n 



Np+i - Np^ 



(3) 



where /i = 0. Some trends with different values for s are represented in Fig. [31 

We want to evaluate the error of the estimated GARCH(1,1) parameters due 
to the difference between the estimated trend and the real one. The statistical 
ensembles for the Monte Carlo simulations are composed by numerically gen- 
erated series composed by a random component and the trend ([3]). First we 
generate a GARCH(1,1) time series {xi} with given parameters {ao, Po, Kq). 
Then we calculate the series ?/„ = Yd=i^iy (corresponding to the logarithm of 
a price series) and we add an automatically generated trend, ^„ = fn + Vn- The 
relation between these two components in the resulting series is characterized 
by the ratio r of the amplitude of the trend and of the noise 

max(/„) - min(/„) 

r = . 

max(?/„) - min(y„) 



If we randomly choose the number s of trend parts between two given values 
■Smin = 1 and Smax = 4 and the ratio r between r^i^ = 0.25 and rmax = 4 we 
obtain a significant statistics. 

From the series {^„} we extract a polynomial estimated trend {/„}, y„ = 
— fn, and we calculate the estimated returns x„ = yn+i — Vn- Then we 
evaluate the GARCH(1,1) parameters {a, (3, K) using the maximum likelihood 
method. When we choose the degree of the polynomial estimated trend we 
must take into account that if it is too large, then the estimated trend begins 
to follow the fluctuations of the noise. From numerical tests it has resulted an 
optimal degree equal to 2s + 3. 
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Fig. 4. The variability of the GARCH(1,1) parameters due to detrending errors with 
respect to the ratio r between the trend and the noise amphtude. The statistical 
ensembles contain 100 generated GARCH(1,1) series with the DJC parameters su- 
perposed on automatically generated trends, (a) The average value of the estimated 
P for different numbers of trend monotonic parts. The continuous line corresponds 
to the value Pdj = 0.8898 used to generate the GARCH(1,1) noise, (b) The rela- 
tive standard deviation of the GARCH(1,1) parameters for the generated noise Xn 
(continuous line) and for the estimated noise x„ (markers) for s = 4. 

4 G ARCH (1,1) parameters variability due to detrending 



Figure H] shows the results obtained by applying the evaluation method of the 
variability of GARCH(1,1) parameters described in the previous section for 
statistical ensembles of 100 time series generated with the DJC parameters 
for different values of r and s. From Fig. it results that the averages of the 
estimated parameter P are randomly distributed around the initial value Pdj 
and they are not influenced by the number of monotonic parts of the trend s 
or by the ratio r. The other G ARCH (1,1) parameters have a similar behavior 
so we have not represented them. Figure Hb confirms this result by means of 
the relative standard deviation for s = 4. Hence, the GARCH(1,1) parameters 
are not influenced by detrending a nonlinear trend if the noise is generated 
using {ttDj, pDj, Kdj). 

In our case the coefficient (3£,j = 0.8898 has a much greater value than a^j = 
0.083. The variability of the GARCH(1,1) parameters at detrending for greater 
ratios q;o//5o is presented in Fig. [51 The number of parts of the artificially 
generated trend is s = 4 and the ratio between trend and noise is r = 2 and 
Q!o and /?o are varied such that ao + Po remains constant, ao + Po = 0.972 
and Kq = Kdj. One observes that for /3o > 0.7 the mean of the estimated 
values i^P^ almost coincides with the initial value Pq and the relative standard 
deviation is less than 3%. Hence, the behavior observed for DJC index is the 
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Fig. 5. The variability of the parameter /? due to the detrending error for statistical 
ensembles of 100 generated GARCH(1,1) series with the DJ parameters superposed 
on automatically generated trends for s = 4, r = 2 and ag + /3o = 0.972. (a) The 
dashed line corresponds to the values f3o used to generate the GARCH(1,1) noise. 

same for smaller values of /?o- But the error (|/?^ — /?o and the relative standard 

deviation a (j3j significantly increases when (3 decreases below 0.6, so in this 
cases the influence of the errors of the estimated trend is significant. The other 
two GARCH(1,1) parameters have a similar behavior. 

Hence the influence of detrending on the variability of G ARCH (1,1) param- 
eters is due especially to the coefficient (3 that generalizes the ARCH model. 
Figure [6] contains the results of a detailed analysis for the minimum value 
= 0.1 in Fig. [5l The variability of the estimated GARCH(1,1) parameters 
significantly increases when the number s of monotonic parts of the trend in- 
creases and the ratio r between the variation amplitude of the trend and the 
noise is larger. 



5 Conclusions 

An important problem in the analysis of financial series is to separate the 
deterministic component and the stochastic one. In this paper we have ana- 
lyzed the influence of the nonstationarity due to a deterministic trend on the 
GARCH(1,1) model and we have shown that for long periods of tens of years 
the results obtained with this model are not sensitive to the existence of a 
nonlinear trend. This behavior occurs for large values of the GARCH param- 
eter f3 and a + f3 close to 1. But for small values of f3 and large values of a 
such that a -|- /? ~ 1 the influence of detrending on the GARCH(1,1) model 
becomes significant. 
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Fig. 6. The mean and relative standard deviation of the GARCH(1,1) parameters 
estimated for x„ obtained on statistical ensembles of 100 time series for each r and 
s. The continuous line in (a) represents ao, in (c) /3o, in (b) represents CT{a) and in 
(d) a{f3). 
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